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Test particles interact with a medium by means of a bimolecular reversible chemical reaction. 
Two species are assumed to be much more numerous so that they are distributed according fixed 
distributions: Maxwellians and Dirac's deltas. Equilibrium and its stability are investigated in the 
first case. For the second case, a system is constructed, in view of an approximate solution. 
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I. INTRODUCTION 

Among the elementary processes of collision between 
ions and molecules, the dissociative attachment of an 
electron to a molecule 

P~ + AB ?± A~ + PB , (1) 

plays a role in the physics of weakly ionized gases 0, • 
However, a full kinetic study is laking. In particular, 
when the electric field is not vanishing, equilibrium solu- 
tions are not available for the distribution functions. 
In the present paper we shall assume that particles 2 and 
4 are much more numerous and can be treated as a neu- 
tral background, while particles 1 and 3 have the same 
(negative) charge e and are subjected to a constant elec- 
tric field E. (Hereinafter, particles P~, AB, A - , and PB 
will be labeled by the subscripts 1, 2, 3, 4. Each of these 
particles is endowed with mass rrii and internal energy of 
chemical link £7j.) 

First of all we recall the full Boltzmann equations for 
particles 1 and 3. The exact form of the collision inte- 
grals is shown. The weak form of the kinetic equations 
is constructed, in order to investigate conservation laws 
and equilibrium. Two cases are considered: 

(1) particles 2 and 4 are distributed according to 
Maxwellians without drift velocity. Under this assump- 
tion, equilibrium and its stability are investigated when 
E = 0. Firstly we show how to construct the equilib- 
rium solutions. Successively, we show the existence of a 
Lyapunov functional for the present problem. A physical 
counterpart of these mathematical results is discussed. 

(2) B is heavier than P _ and A~ . In the limit me — > oo 
particles 2 and 4 can be considered to be distributed ac- 
cording to Dirac's deltas As a consequence, the Boltz- 
mann equations are modified. Moreover, if the electric 
field is small, we introduce a first order spherical har- 
monic expansion for the distribution functions of parti- 
cles 1 and 3. The result is a system of differential-finite 
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difference equations for the problem. By taking advan- 
tage of the smallness of E, a first order solution can be 
constructed for the current functions. 



II. BOLTZMANN EQUATIONS FOR 
PARTICLES 1 AND 3 

The nonlinear integrodifferential Boltzmann equations 
governing the evolution of the distribution function for 
the reacting particles 1 and 3 reads as follows [4| : 

m +v '£ +e %-iL)* = J M +Q M> (2) 

where Ji [f] and Qi [f] are the chemical and elastic collision 

integrals with f = (/i, / 2 , f 3 , U). 

The chemical collision integrals are given by 



JM = J JCi[Qdwdn' 



(3) 



where, for i = 1, we have 

K 1 [§=6(g 2 - m2 )p 3 1 l(g,n-n') 

3 



mi rri2 
raj, m,4 



h{v^)h{w{l)-h(v)f 2 {w) 



,(4) 



being 9{x) the Heaviside step function whereas, for i = 3, 
we have 



34 

mi 7712 



where 



1 

M 
1 

= M 
1 



fMl) h(wll) ~ f 3 (v) ,U(w) 

(mi v + m 2 w + m 4 n') , 

(mi v + m 2 w — m 3 gf% n') , 

12 



'34 = ( m 3 v + m 4 w + m 2 534 re') 



antonio.scarfone@polito.it 



; 34 



M 



(7713 v + 777.4 w — mi g 34 re') 



,(5) 

(6) 
(7) 
(8) 
(9) 
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with M. = m\ + m 2 = rri3 + m 4 . In Eqs. (J3J and (JSJ 
we have introduced the differential collision frequencies 
of the forward and backward reaction ^i|(<?, n ■ n') and 
^Kffi n ' where g = \v — w\ whilst n and n' are 
the unit vectors of the relative velocities before and af- 
ter collision, respectively. We observe that the following 
microreversibility condition holds 

(mi m 2 ) 2 gv\\(g, n-n') = 

(m 3 m 4 ) 2 <7?| ^ 4 2 ( ff , n • n') % 2 - r /34 ) , (10) 

where 
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(9 2 -mj) 



(ii) 



where 



7M$= E ^Cff' n ' n ') \M v le) M w le)-Mv) fi(w) 



and 



1 



«2 = 



(13) 



(rrii v + mi w + mi g n') , (14) 



rrii + mp 

wff = (mjV + miw — mi g n') . (15) 

TOj 4- mi 



with Tjij = 2A4 AE/mi rrij and AE = E 3 + E± — E\ — 
Ei > is the molecular heat of reaction. 

Differently, the elastic collision integrals are given by 



QM = Jn^dwdn' , 



(12) 



The weak form of the kinetic equations for i = 1 and 3 
is obtained by multiplication times a pair of sufficiently 
smooth functions <fii(v) and <j>^v), respectively, integra- 
tion over v, and by summing 



= / £i 



Mv) - MO 



dvdwdn! 



M v lt) } dvdwdn' . (16) 



Observe that for (j>i(v) — 4>z(v) = 1 we get dni/dt + with £ = 2, 4 and [it fixed. 

dn-ijdt = 0, with n, = J" /j dw, that is the total number From the weak form of the Boltzmann equation, by set- 



of test particles is conserved. 



ting 



III. CASE (1) 

We assume particles 2 and 4 much more numerous, so 
that they can be treated as an equilibrium background 
at a fixed temperature T 



h = hi j/i exp 
h = In \ h exp 



> j fi 2 - E 2 - - m 2 v 2 
) { ^ - E± - ~ m 4 v 2 



,(18) 
,(19) 



ft = mf exp 



1 2 

/3 [ - Et - - m e v 



(17) where /j = fi/mf, we obtain 
I 



V = I v\\(g, n ■ n') (mi to 2 ) In 



3l M«S)/4K) r ; 



h{v)h{w) 



h{v) f 2 {w) - / a («S) jkKt)] dudwdn' 



i X]( m * TO ^) 3 / n ■ n') In 



fi(vH)fe(wf e ) r ~ 



fi{v)ft(w) 

r 



[/<(«) /j(t») - //(tog)] dudtodn' < 



(20) 



(i = 1 and 3; £ — 2 and 4), where I? is the left hand side kinetic theory £|, we have 
of (16). Based on these results, by standard methods of 
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Proposition 1. The equilibrium condition 
dh_dh_n 

dt at ' 

is equivalent to 



h{vf t )h{wf t ) = h{v)h{w) , 



(21) 



(22) 
(23) 



with i = 1, 3 and £ = 2, 4 • 
From the first equation we get 



fi = exp 



(3 [ m - ^ - - m, v' 



(24) 



while the second one gives \i\ + [i 2 = M3 + A*4- 

In order to investigate the stability of such equilibrium 
solution, we introduce the following functional: 



C = H + (3 J fx(^x - E x -Jmi« 2 j dv + 13 J f 3 ^ 3 - E 3 - ^m 3 v 2 ^j 

where H = J Tidv, being H = Hi + H 3 and dHi/dfi = ln/j. 

Proposition 2. £ is a Lyapunov functional for the present problem • 
Proof. First of all, from (16) and (20) we verify that 

— = V < . 

dt 

Moreover, by introducing the first Taylor expansion of H around the equilibrium 



dv 



/** = / 



«• + ! f£) * (/. - /■') + (f£) "(A 

H" + a (n - E, - 1 m, i. 2 ") (/, - /,•)/? (/« - £ 3 - J m 3 1. 2 ] (A - /;) 



(25) 



(26) 



(27) 



where * means "at equilibrium", from Eqs. (25) and (27) to Clausius inequality, we observe an additional term due 



we obtain 



dv 



c-c* = J (n-n*y 

+ pj (mi - E X - \ rm v 2 ^j (fx - fx*) dv 
+ P j (ms - #3 - \ m 3 v 2 ^j (f 3 - /*) dv 
= J(H-H)dv . (28) 

Due to the convexity of Ji we can conclude that C > C* . 
The inequality 



dC 
dt 



< 



(29) 



can be interpreted on a physical ground. In fact, by 
introducing the entropy S = —H, we get the following 
thermodynamic inequality: 

dS > — (d£ — dn\ — ^ dn 3 ) , (30) 

where £ = J fx {Ex+mx v 2 /2) dv+ J f 3 (E 3 +m 3 v 2 /2) dv 
is the total energy density of test particles. With respect 



to the fact that the medium of field particles 2 and 4 not 
only provides heat to the gas of test particles 1 and 3 but 
also modifies its composition. 



IV. CASE (2) 

In the limit mg — > oo we have 
mi m 2 mx 
m 3 rrii m 3 
and the following relations hold 



,34 



m 3 



mi 



m 3 



mi 

► w + v~n' , 
w + vn' , 



■ V + Tjx 



w{ 2 — > w 



.34 



w + v + n' 
w — vn' , 
• w , 



where m = 2AE/rrii. Moreover we can pose n 
n' — > $7' and g — > v. 



(31) 



(32) 
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By taking into account that fe(w) = Mi S(w) for £ — 2 and 4, the integrals Ji[£] and Qi[£\ read now 



Ji[f] = J e(v 2 - m )^( v , n-n') 



^i) Ar 4 / 3 («-n')-M/i(«) 

m 3 



d$7 



mi \ 34 , 



m 3 / i! 



i, 12 >+, n-n') 



m3 
mi 



tv 2 n ■ n') + 7V4 n • n')l [/ 3 (« n 7 ) - /*(» 



dn 



(33) 

(34) 

(35) 
(36) 



Equilibrium and its stability for the present problem are 
investigated in §. Our purpose here is to construct 
model equations suitable for an approximate solution. 
As usual in the physics of weakly ionized gases jl| , if both 
the spatial gradients and the electric field are small we 
may resort to a first order spherical harmonic expansion 
of fi(vfl): 



where 



fi(vn) =Ni(v) + n- Ji{v) , 



N i( v ) = 4^/ M vn ) dn > 



(37) 



(38) 



and 



Ji(v) 



47T 



(39) 



By projecting over 1 and CI we get, after some manipu- 
lations, the following system for the new unknowns func- 
tions Fi(0 = Nt(v) and Gift) = J t (v): 



^W 4 F 3 (r)-^ 2 Fi(o 

m 3 J 



-^7- + — V ' G 3(0- — ■^g-^G 3 (0}-{— j v 12(0) (t )W- 



mi 



AA 2 F 1 (e + )-AA 4 F 3 (0 



mi 



mi 



^(D (0 ( ^ ) AA 4 G 3 (r) - A/- 2 Gi(£K 3 2(o) (0 



-7i(0Gi(0 



(40) 



(41) 



(42) 



2e£?_ /? aF 3 (0 _ /'miV ^ 

-7 3 (e)G 3 (o , 



■v^- ... \ 

m 3 d£ \ TO 3 



y 12( 



*> ^ (mr) 3 -^ 1 (^ + ) --^ G 3(o-i 3 2 4 (o) (r 



(43) 



where we have posed ^ = (v ± ) 2 and being f'j^(£) = z^( )(£) — ^im(0 an d 

^k) ® =2w £ ^ v i? & ^ ^ ' ( 45 ) 
with /c = 1 and 2. 

7<(0 = A2 ^(t) (£) + -^4 ^i4(i) (£ ) ' (44) Consider now the stationary space- homogeneous equa- 
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tions We observe that for E = the following equilib- 
rium solutions hold: 



where 



Fi = d exp 



C 3 A/4 / m 3 



'2kTI ' 



exp 



AE\ 



(46) 



(47) 



C1N2 \ m i J ^ry 

(mass action law). Finally, we observe that the electric 



Since are already known, by solving equations i|40|) 
and (|41() . we can obtain the expression of G| , in the 
case of isotropic reaction collision frequency v 



12(0)- 



G« = - 



G« = 



e^Ci V? exp(-mi^/2KT) 

~~k 7 t~ o(t - m )^2 ^ 4 2(0) (0 + 71(6 

e£G 3 



(49) 



'12(0) * 
£ exp(~m 3 C/2kT) 



(50) 




where G^ 1 ' = G^ 1 ' ■ e, with e the unit vector of E. 

In figure 1 we depict, in arbitrary unity, the plots of 



\G ( i\i)\ (full line) and \G\'> {Q\ (dotted line). Since 
the forward equation has a threshold for £ = 7/1, the 
collision frequency of particles 1 suddenly increases, and 
a discontinuity in the relevant plot of |G^ (£)| occurs. 



0-1 1 " • i ■ ■ r " r- 

20 40 60 80 100 

FIG. 1: Plot of the current function {G^ (full line) for the 

charge particles P~ and IG^'^)! (dotted line) for the charge 
particles A - , in arbitrary units. 

field must be a small quantity, \E\ = e, so that we can 
expand Fi and Gi as follows 

Fi = F^+eFjU + ... , Gi = eG [ p + ... . (48) 



V. CONCLUSIONS 



Two linear Boltzmann models have been constructed 
for test particles reacting with a medium of numerous 
field particles. In the first case the field particles are 
distributed according Maxwellians with vanishing drift 
velocity. Theorems on equilibrium and its stability are 
given, as well as their connection with thermodynamics. 
In the second case we consider particles B heavier than 
particles P~ and A - . By means of first order spherical 
harmonic expansion, four equations can be constructed, 
suitable for an approximate solution. 
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